Thermodynamic Limit for the Ising Model on the Cayley Tree 
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While the Ising model on the Cayley tree has no spontaneous magnetization at nonzero tem- 
peratures in the thermodynamic limit, we show that finite systems of astronomical sizes remain 
magnetically ordered in a wide temperature range, if the symmetry is broken by fixing an arbitrary 
single (bulk or surface) spin. We compare the behavior of the finite size magnetization of this model 
with that of the Ising model on both the Sierpinski Gasket, and the one-dimensional linear chain. 
This comparison reveals the analogy of the behavior of the present model with the Sierpinski Gasket 
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The zero-field partition function of the Ising model on 
the Cayley tree is identical to that of the one-dimensional 
linear Ising chain pj , and thus it does not lead to any sin- 
gularities of the zero-field thermodynamic response func- 
tions. This model manifests thermodynamic singular be- 
havior only in response to the magnetic field, and it has 
been well established 0, 0, 0, B that the susceptibility 
diverges in a wide range of temperatures, while sponta- 
neous magnetization remains zero at any finite tempera- 
ture in the thermodynamic limit. 

While the above facts were established three decades 
ago, the exact analytical expression for the zero field mag- 
netization has been derived only recently 0, |(j . Similar 
to the expression for the zero-field partition function, the 
obtained expression for magnetization is equivalent to 
that of the linear chain, except for the fact that the num- 
ber of particles is replaced by the tree generation level. 
It is shown in this work that this fact, together with fix- 
ing an arbitrary single spin, leads to magnetic ordering 
in a wide temperature, range for systems of astronomi- 
cal sizes. This is reminiscent of the effect of extremely 
slow decay of correlations found for the Ising model on 
the Sierpinski gasket [jj , where it has been demonstrated 
that thermodynamic limit represents a poor approxima- 
tion for systems of both laboratory and astronomical di- 
mensions. We compare the finite size behavior of zero 
field magnetization of the Ising model on the Cayley tree 
with that of the Sierpinski gasket and the linear chain. 

For simplicity we consider here explicitly only the Cay- 
ley tree with the branching number (coordination number 



minus one) B = 2, the generalization to arbitrary B be- 
ing straightforward. Following Eggarter Q, we further 
consider a single n-generation branch of a Cayley tree, 
composed of two (n — l)-generation branches connected 
to a single initial site, with the Hamiltonian 
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where J is the coupling constant, H is the external mag- 
netic field, S{ — ±1 is the spin at site i, and (nn) de- 
notes summation over the nearest-neighbor pairs. The 
n-generation branch consists of N n = 2" +1 — 1 spins, the 
O-generation branch being a single spin. The recursion 
relations for the partition function of any two consecutive 
generation branches were found by Eggarter || to be 
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where x = exp(/3J), y = exp(/3if), (3 = 1/ksT is the 
reciprocal of the product of the Boltzmann constant ks 
and temperature T, and and Z~ denote the branch 
partition functions restricted by fixing the initial spin 
(connecting the two (n— l)-generation branches) into the 
{+} and { — } position, respectively. 

It turns out that the nonlinear coupled recursion re- 
lations © cannot be iterated to yield a closed form ex- 
pression in the general nonzero field case Q, 0, 0, 0, • 
However, only recently the present authors have arrived 
6] at the exact analytical expression for the zero-field 
magnetization and susceptibility, by considering the re- 
cursion relations for the field derivatives of the partition 
function, which can be iterated in the limit H — > to 
yield corresponding closed form expressions for systems 
of arbitrary sizes (therefore also for the thermodynamic 
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limit) . The exact expression for restricted magnetization, 
for arbitrary tree generation level n, was found to be 
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where f = tanh(/3J). The same expression was (also 
only recently) independently obtained by Melin et al. [j| , 
using a different approach. 

It turns out that the zero-field magnetization (as well 
as all the other odd field derivatives of the partition func- 
tion) is identically zero in strictly zero field, as the ex- 
pressions corresponding to the two possible orientations 
of the initial spin differ only in sign. However, this is also 
true for an arbitrary Ising system in (strictly) zero field, 
as every spin configuration has a mirror image (obtained 
by flipping all the spins) with exactly the same energy 
and inverted sign of the configurational magnetization, 
and to remove this degeneracy one first has to break the 
symmetry. The usual procedure of breaking the symme- 
try by retaining an infinitesimal field, then taking the 
thermodynamic limit, and finally taking the zero field 
limit, is not suitable here because of the fact that sus- 
ceptibility diverges in a wide temperature region, rather 
than in a single (critical) point. 

Instead, here it seems most suitable to break the sym- 
metry simply by fixing the initial spin in the upward 
position (therefore disregarding all configurations corre- 
sponding to downward initial spin orientation), and in 
what follows we shall first examine the analytical results 
corresponding to this situation. It can be argued, how- 
ever, that the initial spin has a privileged role in this 
lattice, and we shall here also examine numerically the 
consequence of fixing an arbitrary spin (on the surface, 
or in the bulk), on the thermodynamic behavior of ex- 
tremely large (but finite) systems. It will be shown in 
the rest of this paper that fixing any one spin, which is 
equivalent to imposing an infinitely strong local field of 
infinitesimal range (acting on a single spin), leads to sym- 
metry breaking, causing magnetic ordering of systems of 
astronomical sizes in a wide temperature range. 

By fixing the initial spin, one can consider as the order 
parameter the restricted magnetization 
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As t < 1 for T 7^ 0, this expression clearly demonstrates 
that there is no magnetic ordering in the thermodynamic 
limit at any nonzero temperature, confirming earlier re- 
sults 0, 0, IS 3- However, comparing with the corre- 
sponding expression for the one-dimensional Ising model 
(which may also be considered as a Cayley tree with the 
branching number B=l) 
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one sees that the generation level is found in place of 
the number of particles (these two being the same in the 
case of the chain). If one takes the Avogadro's num- 
ber (~ 10 23 ), and the estimated number of hadrons in 
the observable Universe (~ 10 80 ), as the two basic refer- 
ence points in our description of the physical world (as 
opposed to the mathematical definition of infinity), it 
turns out that the slight difference between expressions 
(@J and (J5J leads to strikingly different behavior of the 
two systems. More precisely, it will be shown in this work 
that the Cayley tree remains ordered in a wide temper- 
ature region for systems by far exceeding, in number of 
particles, the size of the observable Universe. This be- 
havior is reminiscent of that of the Ising model on the 
Sierpinski Gasket, and may provide some further insight 
into the thermodynamics of model systems on hierarchi- 
cal and fractal lattices in general, the latter nowadays 
being commonly accepted as frequent representatives of 
our physical world. 

Let us begin this (extreme) finite size analysis by ex- 
plicitly plotting magnetization of the Cayley tree and 
the Ising linear chain, as given by expressions J2J and 
©, for various system sizes N < 2 2 "' +1 — 1 ~ 10 309 
(N = 2 n+1 - 1 for the Cayley tree, and N = n+l for the 
Ising linear chain), as shown in Figs. la) and lb), respec- 
tively. It may be interesting to mention here that the 
extremely simple form of expression |J5J is deceptive in 
the sense that it requires special numerical handling when 
such large system sizes are considered (for N = 2 2 +1 — 1 
two thousand digit precision was needed by the Maple VI 
software for algebraic manipulation, for numerical calcu- 
lation of magnetization at low temperatures). 

One can see from Fig. la that even for these incon- 
ceivably massive systems (with more then 200 orders of 
magnitude more particles than the estimated number of 
hadrons in the observable Universe), magnetization of 
the Cayley tree remains substantially removed from its 
limiting value of zero, while from Fig. lb it is seen that 
the magnetization of the Ising linear chain for the same 
system sizes far more rapidly collapses to zero, and can- 
not be distinguished from its limiting zero value on the 
scale of the graph. 

We have also performed similar calculations for the 
Ising model on the Sierpinski Gasket, for which it was 
established by Liu Q that the decay of correlations is ex- 
tremely slow, and the mean square magnetization (m 2 ) n 
approaches its limiting zero value in such a slow fashion 
that the thermodynamic limit represents a poor approx- 
imation for both laboratory size systems, and systems of 
astronomic dimensions. In this case the lattice at stage n 
is constructed by assembling three triangular structures 
of stage n — 1, and rather then considering ^m 2 ) as was 
done in reference Q, one can break the symmetry by 
fixing the three vortex spins in the current construction 
stage, and consider the corresponding restricted magne- 
tization < m >+ ++ . In Fig. 2 we show the comparison 
of the points T m , where the restricted magnetization falls 
to half of its maximum value, for the Ising model on the 
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FIG. 1: Magnetization of the a) Ising model on the Cayley 
tree, and the b) the Ising linear chain, for systems with N = 

2 2 +1 — 1 particles, with k = 1,2, ...,10. Open circles represent 
the points where magnetization falls to half of it's maximum 
value. 



Sierpinski gasket, the Cayley tree, and the linear chain. 

It is seen that the magnetization of the Cayley tree 
demonstrates behavior similar to the Sierpinski gasket, 
with an extremely slow decay of ordering with the num- 
ber of particles (observe the double logarithmic scale on 
the abscissa), in complete contrast with the linear Ising 
chain. It should be noted that while both in the case of 
the linear chain and the Sierpinski gasket the considered 
system sizes also imply astronomic linear dimensions, in 
the case of the Cayley tree the system linear dimension 
remains miniscule (for the largest considered systems of 
N = 2 2 +1 — 1 particles, the largest distance between any 
two spins remains only 2 10 = 1024 lattice units). In the 
case of the Cayley tree one cannot therefore talk of "long 
range correlations" , rather, the extremely slow decay of 
magnetic ordering with temperature may be attributed 
to the short system span and the fact that every spin 
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FIG. 2: Points where magnetization falls to half of its maxi- 
mum value, for the Ising system on the Sierpinski gasket (full 
squares), Cayley tree (open circles) and the linear chain (full 
circles), as a function of the number of particles. The two ver- 
tical dotted lines indicate system sizes corresponding to the 
Avogadro's number, and the estimated number of hadrons in 
the observable Universe, respectively. 



"sees" a huge number of other spins a short distance away 
(perhaps a term "wide range correlations" would be more 
appropriate) . In the case of other fractal and hierarchical 
lattices, it remains to be seen to what extent the relation- 
ship between "long range" and "wide range" correlations 
may contribute to the magnetic ordering of finite sys- 
tems. In particular, the question remains whether real 
porous samples with magnetic properties that may well 
be microscopically explained by model systems which do 
not display spontaneous magnetization in the thermody- 
namic limit, still may exhibit magnetic ordering at mea- 
surable temperatures, due to both their intricate struc- 
ture and their finite size. 

As mentioned earlier, it can be argued that the initial 
spin occupies a privileged role in the Cayley tree, and 
here we address the question of whether fixing any other 
spin is sufficient to break the symmetry, and cause order- 
ing of massive systems in a wide temperature range. Let 
us therefore consider the situation of building the Cayley 
tree in the regular fashion (using recursion relations (J2J)) 
up to level p, fixing a single spin at that level, and then 
applying modified recursion relations (corresponding to 
asymmetric restricted partition functions) up to level q. 
From here on, we shall term the branch containing the 
fixed spin as the "left" branch, denoting the correspond- 
ing restricted partition functions by w„ , in parallel with 
the regular "right" branches (that do not contain any 
fixed spins), with the "regular" restricted partition func- 
tions . Choosing p — corresponds to a single surface 
spin being fixed, p = q corresponds to fixing the initial 
spin, and < p < q corresponds to fixing a single bulk 
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spin at level p. 

Formally differentiating equation (J2J with respect to 
field, the recursion relations for the field derivatives of 
the restricted partition functions up to level p are found 
to be @ 
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where we have used notation ft = f3H for the reduced 
field. 

At all following levels (including p+ 1) the "right" par- 
tition functions will differ from the "left" partition func- 
tions, and the general recursion relation is given by 

wt +1 = y ±x [K^ 1 + w-x* 1 ) {Z+x* 1 + Z-x^)] , 

(7) 

where for the derivatives we have 
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We now take one "left" branch of level p (by setting 
w+ = and Wp = 0), and one "right" branch of level 
p (using both Z+ and Z~), to form the "left" branch 
structures at level p+ 1 (if p = q we do not continue with 
forming the p + 1 level structure, rather just retain Z+ 
as the final partition function). The recursion relations 
@ and ijSJ) for the restricted partition functions an their 
derivatives, respectively, at this step take explicit form 
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In the following steps (p + 1 < n < q), the "left" 
branches are combined with (regular) "right" branches 
to form new "left" branches through general recursion 
relations Q and |(HJl, while the "right" branches are all 
the time calculated using J2J and iJBJ. 

Evidently, the set of equations J2J ? 1|6I1U[1 is too com- 
plex to be iterated analytically into a closed form ex- 
pression, for arbitrary p and q. On the other hand, using 
algebraic manipulation software such as Maple or Math- 
ematica, it is a simple exercise to perform these iter- 
ations numerically for arbitrary size and temperature, 
with arbitrary precision (depending only on the avail- 
able computer speed and memory, where present day 
personal computers are more then well suited for the 
job). We have performed such calculations for various 
values of p, for q < 1024 (corresponding to systems of 



up to 2 2 +1 — 1 ~ 10 309 particles), in all cases it was 
found that fixing any spin indeed does break the sym- 
metry, inducing ordering of enormous systems in a wide 
temperature range, while fixing the initial spin does pro- 
duce a somewhat stronger effect then fixing a single sur- 
face or bulk spin. In fact, it is found that temperature 
T m , where magnetization falls off to half of its maximum 
value, monotonously increases with p, for any given gen- 
eration level q. This effect may be attributed to the fact 
that looking from a given position, the average distance 
to all the other spins monotonously decreases as the po- 
sition shifts from surface towards the center. Larger dis- 
tance from the point of fixing a spin implies a weaker ef- 
fect, in accordance with the numerically observed behav- 
ior. In Fig. 3 we show the comparison of the T m curves 
as a function of generation level q for p = (full circles), 
corresponding to fixing a single surface spin, and p = q 
(open circles), corresponding to fixing the initial spin. 
From the figure, it follows that the difference between the 
values of T m obtained by fixing the initial spin and by fix- 
ing a single surface spin diminishes with system size, and 
that the conclusions about ordering of extremely large 
systems in a wide temperature range continues to hold 
(as well as the comparison of the current model with the 
Sierpinski gasket and the chain), if one chooses to break 
the symmetry by fixing an arbitrary spin. 

In conclusion, as the Cayley tree in itself (with its infi- 
nite dimension and nonzero fraction of surface nodes) 
represents a highly unphysical system in the thermo- 
dynamic limit, modeling real world physical phenom- 
ena (such as, for instance, blood vessels or the nervous 
system) requires more attention focus on the finite size 
behavior. In this work it is shown that the thermody- 
namic limit represents a questionable approximation for 
the Ising model on the Cayley tree, as the magnetic or- 
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dcring decays extremely slowly with the number of parti- 
cles. As the susceptibility diverges in a wide temperature 
range in the thermodynamic limit, and here we study fi- 
nite (albeit extremely large) systems, symmetry breaking 
is achieved by fixing a single (any one) spin, rather then 
by the usual procedure of retaining an infinitesimal field 
while taking the thermodynamic limit, and only then tak- 
ing the limit of zero field. On a double logarithmic scale 
with respect to the number of particles, we have found 
that Ising model on the Cayley tree behaves in a fashion 
rather similar the behavior of the Ising model on the Sier- 
pinski gasket, in contrast with the Ising chain. However, 
as opposed to extremely long range correlations that are 
present in case of the gasket, this behavior may be at- 
tributed to extremely "wide range" correlations. 



FIG. 3: Temperature T m where magnetization falls to half 
of its maximum value, for the Ising system on the Cayley 
tree, as a function of system size N = 2 q — 1, on the double 
logarithmic scale. The upper curve (open circles) corresponds 
to fixing the initial spin (p = q), and the bottom curve (full 
circles) to fixing a single surface spin (p = 0). For all the 
other values of < p < q, T m is found to lie between the two 
shown curves, monotonously increasing with p. 
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